Dark matter seeding and the kinematics and rotation of neutron stars 
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Self-annihilation of dark matter accreted from the galactic halo in the inner regions of neutron 
stars may affect their kinematical properties, namely velocity kicks and rotation patterns. We find 
that if a stable long-lived single or multiple strangelet off-center seed forms, there is an associated 
change in momentum and torque that may affect the kinematical observables of the star. 



Light DM particle candidates in the range of 4.5 — 12 
GeV/c 2 are favoured by recent results in DM direct 
detection, including the DAM A/LIBRA, CoGENT and 
CRESST-2 experiments [1], although these results have 
yet to be confirmed. Some DM indirect detection exper- 
iments Q support these indications. There is however 
tension with other results, notably Fermi dwarf limits 
3 and CMB damping tail measurements but these 
can be avoided if for example there is an asymmetric 
DM component [5J. In addition, these particles interact 
gravitationally and can be accreted onto stars. Here we 
explore some implications for neutron stars (NSs) that 
extend our earlier work [6]. 

By studying the physics of NSs, we find that macro- 
scopic measurable quantities such as linear velocities and 
rotation periods could provide an indirect indication of 
the presence and nature of DM. Typically, the average 
spatial velocity of a NS progenitor is less than 15 kms -1 , 
which means that they must receive a substantial kick 
at birth. Pulsars (PSRs) show very regular electromag- 
netic emission due to the misalignement of rotation and 
magnetic axes that make them visible to the intersect- 
ing light-of-sight observer on Earth. Surface magnetic 
fields for these objects are believed to be in the range 
B ps 10 9 - 10 15 G. As is well known, most radio PSRs 
exhibit transverse velocities ranging from a few 100 km 
s _1 to above 1000 km s -1 , mostly in the galactic plane 
Q. As for rotation, the observed PSR periods range from 
less than 0.1 s to the death line at « 1 s. A different origin 
accounts for the recycled millisecond PSRs or magnetars 
with somewhat larger periods [8[. It has been shown in 
recent simulations [9J that rotation and kicks at birth can 
be linked to the detailed mechanisms in supernova (SN) 
core collapse and this, unless the kick is head-on, may 
cause a change in the rotation pattern. One can obtain 
modest kicks in the 100 — 500 km s~ 4 range [10] based 
on an anisotropic momentum distribution of the ejecta 
or parity- violating scattering of neutrinos [Tlj |. 

Previously @ , we showed that the effect of DM seed- 
ing in NSs may provide a mechanism to form strange 



quark matter (SQM) bubbles or strangelets that, if sta- 
ble and long-lived, could induce a partial NS conversion 
into a hybrid SQM star. This effect is somewhat analo- 
gous to what happens in nuclear ignition models where 
sparks can be induced in the bulk by heating a nuclear 
plasma with a laser [13]. Terrestrial strangelet produc- 
tion searches have been performed at high energy collid- 
ers at E864 in BNL and are currently undertaken at LHC 
and ISS (International Space Station) 13j. We show 



here that DM single or multispot spark seeding in com- 
pact objects may lead to a modification of NS kinematics 
and rotation. A typical NS has a mass Mms ~ 1.4M Q , 
radius Rns ~ 10km and velocity, vns, that can be 
in excess of ps 1000 km/s. The associated momentum 
p NS w M NS v NS ps 3 10 33 g 10 s cm/s ps 3 10 41 g cm/s. 
The birth energy released in a SN type II event is known 
to be E v pa 10 53 erg. If kicks were to be explained by 
(massless) neutrinos, this means that they should carry 
typically p„ ps E v /c ps 10 53 /3 1 10 = 3.3 1 42 g cm/s. 
This results in an asymmetry a p = pns/Vu ~ 10 2 ■ 
However, this neutrino rocket mechanism may be resid- 
ual since thermal effects tend to make emission isotropic. 
Recent analyses of PSR observations seem to indicate 
possible correlations for the directions of the spin axis 
and the velocity vector, another piece of evidence for a 
large velocity in 3D phase space [l4| . The radial distribu- 
tion of the birth location of NSs at given galactocentric 
distance R can be fit using a relation derived by Lorimcr 
that includes selection effects UM 
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where n = 2.35, R = 1.528 kpc and A = 64.6 kpc" 4 35 
and peaks at R ma x = 3.6 kpc. NSs can accrete DM 
from the galactic halo by gravitational capture. This 
DM profile can be characterized using the continuously- 
varying slope function |16j |. 



Pdm{t) = p- 2 e 



(2) 



with p-2 = 0.22 GeV/cm 3 , a — 0.19, and r_2 = 16 kpc so 
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that at the solar neighbourhood the Keplerian velocity is 
v w 220 km/s and the local DM density is 0.7GeV/cm 3 . 
Both observables can be approximately fit assuming most 
of the matter is DM and using typical parameters of 
Milky way-size halos. DM is accumulated over time in 
the central NS regions where it is rapidly thermalized 
(t t h ~ 0.13 ps to achieve T — 30 MeV or t t h ~ 4/is to 
achieve T = 1 MeV 17] in the hot or cooled NS, respec- 
tively). Assuming DM is well described by a Boltzmann 
distribution, the thermal radius is given by 



r«,(t) 



3k B T c (t) 
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where nix is the DM particle mass. This quantity is 
time-dependent since the central temperatures T c and 
mass densities p c evolve with those in the progenitor 
star. For example, during the pre-collapse, typically 
T c = 0.001 - lMeV and p c = lO 10 " 13 g/cm 3 while in 
the bounce phase, the temperature increases to T m 40 
MeV and there is a dramatic increase in the central den- 
sity up to p c — 4 10 14 g/cm 3 in the time scale of the 
gravitational collapse, t co u apse 10 -3 s. Assuming DM 
is of the Majorana type, it may self-annihilate. We de- 
fine a r = n th as a measure of the radial coordinate 
fraction involved in annihilation processes that, at the 
few GeV scale, involve photons, leptons , light qq pairs 
[l8j |. The ratio of the thermalization volume to that in 
the NS, Vth/Vffs ~ a 3 w 10~ 6 . The central annihila- 
tion or spark production rate at radius r will be given 
by dr J^ r ' ) ps (<J a v) n x 4:irr 2 where we assume that the 
DM particle number density, nx(r) ~ nx, is constant 
on a scale r « 10 2 m. In a SN type II core collapse, 
peak temperatures and densities are achieved off-center 
and production of one or more sparks may subsequently 
happen in those regions. Other exotic mechanisms may 
induce thermal or quantum nucleation of quark bubbles 
[lij . however, it turns out that this may not be efficient 
since typical timescales are larger than standard cooling 
scenarios. The energy released over the thermalization 
volume and, accordingly, the rate of spark formation can 
be obtained assuming an efficiency, /, for DM annihila- 
tion @ as E a « 2fCm x c 2 where C is the DM capture 
rate. Assuming NS matter is mostly in the core under 
deconfined nucleons, the capture rate can be written as 
201. 



2.7 x 10 29 p DM 
m x {GeV) po 



(4) 



where we have used cjxn = 7 10~ 41 cm 2 and average core 
density pns = y-* s ■ We consider the spatial NS dis- 
tribution speak and the Keplerian velocity is calculated 
from the enclosed DM and baryonic (B) mass compo- 



nents as v{R max ) 



_ . I GJADM + BiRn 
Rmax 



Galaxy mod- 



Mb/Mdm ~ 1 , then at location r, MoM+B(r) ~ 
Mdm{t) + Mb = Jo ^r' 2 p DM (r')dr' + M B where M B = 
J Pb^ This yields v 220 km/s at the solar circle 
and at the NS distribution peak v ps 190 km/s. The 
spark energy release AE spar k = 2/to^c 2 will be de- 
posited in the medium in the vicinity of the production 
site. The primordial spark energy ratio to nucleon mass 
(tojv f» lGeV) is AE s /m N « 2m x f/m N w 0(1 - 10). 
This high energy deposition may potentially lead to the 
formation of a strangelet @ . The binding energy of such 
a strangelet is given in the MIT model as [22| E^ ct « 
E A (pi,mi, B) + Ecoui where pi and m s are the chemical 
potential and mass of the ith-type quark, respectively. B 
is the MIT bag constant and Ecoui is the correction due 
to electrical charge. Those with mass number A > A m i n 
are long-lived and energetically stable. As an estimate, 
the production rate of long-lived strangelets (with A = 
10) at R max would be 7V slct = E a /E^ ct « 10 27 s' 1 for a 
density of about 2no,A (no, a = 0.17 fm -3 or equivalently 
mass density po — 2 10 14 g/cm 3 ). The ratio of spark en- 
ergy to strangelet binding energy BE — E^ ct — ArriN is 
AE S /BE w 0(1 — 10) making strangelets energetically 
accessible systems. 



els [21| show an inner mass distribution where / 
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We will consider now for simplicity a single spark seed 
model located off-center at r ff = r //e ff, r ff < 
ct r RNS- From the typical conditions in the interior 
of the (proto-)NS in the location of the growing spark 
a r r; 10~ 2 — 10 _1 . The spark radial direction forms 
an angle /3 with respect to e sp i n , unitary vector in the 
spin axis direction. Once formed, the quark bubble is 
composed of burnt material consisting of u, d and s 
light quark sectors. The seed radius R s will grow as a 
function of time, t, and central distance, r, due to the 
strong and weak interactions at a rate dR s /dt — Vb(t,r) 
that can be estimated assuming a simplified two — step 
burning front that could evolve dynamically 23] depend- 
ing on the matter equation of state (EOS). Since the 
strong interaction time-scale t str0 ng < C 1 , the annihi- 
lation mechanism may produce an effective local heat 
source. Assuming the minimum stable strangelet has 
baryon number A m i n , then in a semiclassical approach 
there is an increase in the energy release per nucleon, 
Ae th w £^ « « 0.1 - lGeV w 0(A QCD ) 

which is of the order of the binding energy of the nucleon 
or QCD scale, A QCD ~ 0.217 GeV. 

In order to obtain the effect on the kinematic NS 
observables, we consider the combustion of an electri- 
cally charge neutral system of hadronic unburnt matter 
(phase 1) into the burnt quark phase (phase 2). For the 
hadronic star (composed of neutrons, protons and elec- 
trons in beta equilibrium and in global electrical charge 
neutrality), the EOS P = P(e,p) will be taken as the 
APR model (ill at T 0, No magnetic fields will be 
considered. Instead, for the deconfined quark burnt ma- 
terial the MIT bag model will be taken for consistency 
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with the finite lump of quark matter description. As the 
combustion takes place, weak interactions act on a time- 
scale of t wea k ~ 1 /r Kaon ~ I0~ 8 s and a u — d — s — e 
system forms while the gravitational force tends to con- 
tract the star on a timescale t grav « 2i?jvs/c = I0~ 4 s. 
Non-zero strangeness fraction involves reactions includ- 
ing u + d — > u + s, it + e — s> e? + i/ e , u + e — > s + u e , 
where neutrinos freely escape. Jump conditions near the 
front in this relativistic fluid are given by a set of con- 
servation laws (energy-momentum and baryon current 
) across the combustion front and the increase in en- 
tropy condition for the two phases. In the rest frame 
for the burning front (c = 1) [25[ the set is written as 
(ei+pi)v'l+pi = (e 2 +p 2 )vj+p 2 , (ei+pi)v! = (e 2 +p 2 )v 2 , 
siVi < s 2 v 2 , mvi = n 2 v 2 , where Pi,ei,m = Pi/mi and 
Vi are the pressure, energy, number density and veloc- 
ity for the ith-phase. Using the preceding equations 



results in 



,,2 _ (P2-Pl)(e2+Pl) 



and 



,,2 _ (P2-Pi)(ei+P2) 



1 («2— «i)(«i+P2) 2 (c2— ei)(e 2 +pi) ' 

It is possible to classify the conversion when v\ , v 2 are 
comparable to the sound velocities for the two phases, 
c li = if 1 - According to accurate calculations [26|], fi- 
nite strangelets may be positively charged on their sur- 
face due to the strangeness content in their deep interior 
for A « 10 2 — 10 s and cannnot bridge the gap to large 
A systems. This may indicate a limited growth poten- 
tial in the context of safe finite volume ultrarclativistic 
heavy ion collisions in vacuum, however, in the core of 
these astrophysically-sized scenarios, the pressure is very 
large and hybrid objects may exist. The radial veloc- 
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FIG. 1. Hadronic, v\, and quark matter, V2, burning velocities 
as a function of matter mass density. The burning front has 
a velocity vj, = —vi. 



ity profile in the star can be constructed by solving the 
static Tolman equations [27j or their rotating counter- 
parts. We will not consider the more realistic rotating 
profile and, instead, we will use the r-dependent velocity 
in the burning site as in 23[ as v(r) — jpdr' . Once 
formed, the growth of the quark lump can build up a 
net excess of momentum AP in the system as a result 
of the different pressures in the burnt and unburnt re- 



gions, AP12 = Pi — P 2 - The associated impulse, Z, and 
force, AF relate as I = AFAt = AP with At the time 
over which the thrust is acting. We make the simplify- 
ing assumption that the burning front remains spherical 
and moves outwards at Vb(r) = — i>i(r) and reaches the 
edge R at At = f R -^f^ in the comoving frame. Due 

° Jr off v b (r) o 

to the fact that the lump is evolving in a non-inertial 
frame there will be associated forces and, for simplic- 
ity, we adopt the weak field Newtonian approximation. 
The mechanical force acting in the radial direction can 
now be parametrized using the effective transverse area 
AS = nR 2 s {t) as AF mech = AP 12 ASe off . There will 
be a gravitational acceleration component due to the un- 
burnt mass M(r) acting on the quark lump with vol- 
ume V s — t^itRs and mass M S (R S ) = V s (R s )p 2 (r f f) and 



therefore AF 



.„■„,. - M s (i? s )g where g = - GM Jf off) e o// , 



using the fact that the unburnt progenitor star compo- 
nent is M(r ff) = V(r ff)pi(r ft). Due to the fact the 
compact star is rotating with angular velocity 1~2, there 
is an additional effect due to centrifugal force AF centr if 
as the quark lump grows towards the star edge. This 
can be estimated as AF cen t r if = —M S (R S ) f! x x r ff . 
Finally, AF centr i} = V S (R S ) p 2 (r ff) fl 2 r fj Sin/3. Ad- 
ditionally there will be a contribution from the Corio- 
lis force AF cor due to the motion of the burning front 
that will induce some extra rotation. The form adopted 
is AF cor = -2M s (R s )Sl x v b . That is AF cor = 
V S (R S ) p 2 (r ff)Qvi(r ff). For the sake of simplicity we 
will now consider motion in the equatorial plane and 
set j3 = 7r/2. Then, the composition AF = AF mec h + 
AF grav + AF centr i / will produce an impulse that can im- 
part a velocity kick in the progenitor NS with mass M^s 
as, 



Av NS = 



X 



M 



NS 



AP 

M 



NS 



(5) 



Then the combustion will partially (or fully) exhaust the 
remaining unburnt material and then AF — > 0. As- 
suming that in the final configuration, the object has 
a moment of inertia If the angular velocity is affected in 
the conversion, since there must be a torque r modifying 
the angular momentum A J = rAt = r Q ff x AP. Using 
AF = AF cor these expressions we finally derive, 



Afl 

IT 



ArAt 



(6) 



The moment of inertia If will be, in general, time- 
dependent and smaller than the original due to EOS soft- 
ening in the conversion. As an approximate parametriza- 
tion to this value, we take If — dMNsR^s with a 
d ~ 0.18. It is important to note that if the lump is not 
able to grow within the star structure due to energetics 
then the star configuration could be hybrid 'J^S • This 
event could provide an efficient source of gravitational 
waves and set a wave-form pattern with clear signatures 



4 



of the conversion. The possible contribution to the jet- 
like formation in the initial stages of a NS could be some- 
what influenced by internal shocks and crust ejection if 
the conversion is complete [29[ . 
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FIG. 2. Average velocity kick in km/s (upper panel) and 
relative change in angular velocity (lower panel) as a function 
of hadronic mass density, p. 



In Fig. [T]we can see the allowed physical velocities of 
the two phases as obtained for matter in beta equilib- 
rium using the hadronic APR EOS [HI and the MIT bag 
model EOS for the quark phase as a function of the loga- 
rithm of the mass density. We use an intermediate value 
for the bag constant range, B 1 / 4 = 150 MeV, and T = 0, 
where long-lived strangelets can form. The solutions are 
constrained to the physical values < v\ < 1. It can be 
seen that for densities logi [p(g/cm 3 )] < 14.685 burning 
is not allowed. The particular density window could be 
EOS-dependent but we expect this feature to be a gen- 
eral trend. In Fig. [2]the variations in linear velocities Av 
in km/s (upper panel, dashed line) and angular velocity 
ratio Afl/fl (lower pannel, solid line) are shown as a func- 
tion of mass density. We have assumed a typical initial 
value for frequency of v = 300 Hz. We see that the linear 
velocities obtained are indeed bimodal, in the sense that 
both high kicks and low kicks are possible and are cor- 
related with the change in the angular velocity. These 
kicks are compatible with NS observed values and also 
with measurements of the change in the angular veloc- 
ity during GRB emission where Newtonian estimations 
give Afl/fl w 10~ 3 - 10~ 2 [3(J. We have used on aver- 
age tthrust ~ 0.3 ms according to radial profile estimates 
of r; 1 km burning propagation in the allowed density 
window. We have shown that if a stable quark seed is 
able to grow, due to the annihilation of DM particles in 
the central regions of compact objects and according to 
the EOS of matter, then a kinematical signature should 
be present giving both bimodal low and high kicks and 
spin angular velocity values. We would like to thank J. 
Miller and C. Albertus for helpful discussions. M. A. 
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